Environment-dependent dissipation in quantum Brownian motion 
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The dissipative dynamics of a quantum Brownian particle is studied for different types of environ- 
ment. We derive analytic results for the time evolution of the mean energy of the system for Ohmic, 
sub-Ohmic and super-Ohmic environments, without performing the Markovian approximation. Our 
results allow to establish a direct link between the form of the environmental spectrum and the 
thermalization dynamics. This in turn leads to a natural explanation of the microscopic physical 
processes ruling the system time evolution both in the short-time non-Markovian region and in the 
long-time Markovian one. Our comparative study of thermalization for different environments sheds 
light on the physical contexts in which non-Markovian dissipation effects are dominant. 

PACS numbers: 03.65.Yz, 03.65.Ta 



I. INTRODUCTION 

The study of open quantum systems has recently 
received renewed attention due to the importance of 
environment-induced effects, such as quantum decoher- 
ence, both in fundamentals of quantum theory and in 
newly emerging quantum technologies [l|, 0] ■ Quantum 
systems are never completely isolated from their external 
environment. The interaction between the system and its 
surroundings induces decoherence phenomena destroying 
quantum superposition and entanglement. Environment 
induced decoherence has been studied extensively for the 
damped harmonic oscillator model. In this context, the 
decoherence of a superposition of Gaussian wave packets, 
a prototype of Schrodinger cat state, has been studied 
theoretically and experimentally 0-0] ■ 

Decoherence also plays a major role in quantum in- 
formation technologies since the coherence time typically 
determines the operational time of a quantum device, 
e.g., a quantum logic gate. Understanding the dynam- 
ics of exemplary quantum systems interacting with their 
environment is therefore of crucial importance for both 
fundamentals and applicative aspects of quantum theory. 

In this paper we investigate the dynamics of a damped 
harmonic oscillator, or harmonic quantum Brownian mo- 
tion (QBM) model, interacting with different thermal 
bosonic environments. The QBM model is one of the 
few models of open quantum systems amenable to an 
analytical solution [11,0, 

A very general derivation of the nonlinear Langevin 
equations for a damped harmonic oscillator, for general 
microscopic system-environment couplings, is given in 
Ref. [2^. In the case of bilinear system-environment 
coupling the generalized exact master equation for the 
reduced system is known as the Hu-Paz-Zhang master 
equation [l3l|. This master equation is typically solved 
numerically. In some cases, explicit solutions in closed 
form exist, e.g., for an initial Gaussian wave packet or a 
superposition of Gaussian wave packets [1^. In this pa- 
per we study the dynamics for initial Fock states of the 
harmonic oscillator and we use a perturbative approach 
that nonetheless allows to study non-Markovian features 



due to structured environments. Our aim is indeed to 
obtain simple analytical expressions in closed form for 
the observables of interest in order to gain insight in 
the fundamental microscopic processes ruling the non- 
Markovian dissipative dynamics. 

The QBM model is widely used in many physical con- 
texts. Indeed it describes a quantum electromagnetic 
field propagating in a linear dielectric medium '27], a 
particle interacting with a quantum field in dipole ap- 
proximation pL9| and a single trapped ion subjected to 
artificial colored noise [l^ . In addition to these quantum 
optical applications, the QBM model is used in nuclear 
physics [_28j and quantum chemistry [2^. For this reason 
the literature on this model is vast and crosses several 
fields of science. Interesting results have recently led to a 
better understanding of anomalous diffusion for the free 
QBM model, i.e., in absence of a trapping potential [3Q| . 
The dynamics of an initial Gaussian state in an anhar- 
monic potential have also been studied [sij . 

The Hu-Paz-Zhang master equation has recently been 
used to study the entanglement dynamics of initial co- 
herent and twin-beam states of two non-interacting har- 
monic oscillators linearly coupled to common |32h34| | or 
independent structured reservoirs [ssl l36j . For inde- 
pendent reservoirs the entanglement dynamics for initial 
non-Gaussian states was presented in [s^- The non- 
Markovian dynamics of entanglement for two coupled 
harmonic oscillators was investigated in for differ- 
ent types of environment. Very recently strategies of op- 
timal decoherence control have been demonstrated for 
non-Markovian two-level systems [3.9J. We focus in this 
paper on a system more complicated than a two-state 
system, i.e., a single quantum harmonic oscillator, but it 
would certainly be of interest to extend the optimal con- 
trol analysis to systems with a non-finite Hilbert space. 

In this paper we use a time-convolutionless perturba- 
tive master equation that does not rely on the Markovian 
approximation and can therefore describe situations in 
which the spectrum of the environment has a structure. 
This is, e.g., the case of atom lasers 14011 or atoms decay- 
ing in photonic band gap materials [4111. By specifying 
the form of the spectrum we obtain analytical expres- 



2 



sions for the mean energy of the system in a close form. 
In this way we can establish a clear connection between 
the reservoir spectral properties and the non-Markovian 
dynamics of the quantum Brownian particle. 

The coupling between the system and the quantized os- 
cillators constituting the bosonic reservoir is given, in the 
continuum limit, by the reservoir spectral density. Differ- 
ent physical contexts are characterized by different forms 
of the reservoir spectrum. The three main classes typi- 
cally considered in the literature are the so-called Ohmic, 
sub-Ohmic and super-Ohmic spectra. We will compare 
the thermalization process for these three types of reser- 
voirs. 

It is worth mentioning that recent advances in reservoir 
engineering techniques [6] pave the way to experiments 
aimed at simulating paradigmatic models of open quan- 
tum systems as the one considered in this paper. In the 
trapped ion context, e.g., the simulation of a QBM model 
for an Ohmic environment has been proposed in Refs. 
(2^ [42} . The same method can be extended straight- 
forwardly to simulate the sub-Ohmic and super-Ohmic 
environments here considered. These experiments would 
allow to test in a controlled way a fundamental and ubiq- 
uitous model such as QBM. Understanding which type 
of environment leads to the faster or slower decoher- 
ence/dissipation dynamics can be of great importance 
in the choice of the physical system for implementing re- 
alistic quantum devices such as a quantum computer. 

The paper is structured in the following way. In Sec. 
we introduce the model under study and the master equa- 
tion describing the dynamics. Section Hill introduces the 
three different examplary reservoirs, namely the Ohmic, 
sub-Ohmic and super-Ohmic reservoirs, used in our com- 
parative study. In Sees. IIVI and IVl we discuss our results 
for the decay rates and heating dynamics of the QBM 
model. Finally, in Sec. IVII we present conclusions and 
outline possible future prospects. 



II. MASTER EQUATION FOR QBM 

Let us consider a quantum particle of mass m moving 
in a harmonic potential. The Hamiltonian of the system 
is 



Hs = loq a^a + 



1 



(1) 



where a and are the creation and annihilation oper- 
ators of the quantum harmonic oscillator, loq is the fre- 
quency and h is set to 1. The environment is a heat bath 
modeled as an infinite chain of harmonic oscillators 



n=0 ^ ^ 



(2) 



where 5„ and b]^ are the creation and annihilation opera- 
tors, respectively, and a;„ is the frequency of the nth oscil- 
lator. The system and the reservoir are coupled linearly 



via the position operators, X oc a + and Xn oc b + 
for the system and reservoir oscillators, respectively, so 
that the interaction Hamiltonian is given by 



Hi = -4=(a-}-a'l')y^fc„(&„ 



bV 



(3) 



where fc„ measures the coupling between each reservoir 
mode and the system oscillator. 

A master equation describing the QBM dynamics can 
be derived starting from the total Hamiltonian 



H = Hs 



aHr 



(4) 



where a is a dimensionless constant proportional to the 
strenght of the coupling between the system and the en- 
vironment. In the weak couling limit (i.e., when a 1), 
assuming initially factorized state {p = ps ® Pe) and a 
thermal reservoir, we obtain the following secularly ap- 
proximated master equation for the damped harmonic 
oscillator [l| 



dt 



Ps{t) 



A(t) - -fU) . t . 
\2a^ psa - aa'ps - Psaa' j 



A(t)+7W 



{2apsa^ - a)aps - Pso^a) , 



where 

m 



-2 I dt' I dujj{uj) 
Jo Jo 

X cos{ujt') cos{ujQt'), 



Niio) + 



(5) 



(6) 



7(t) =2 / dt' 

Jo Jo 



> tj(uj') 

duj — - — sin(cL't') sin(a;oi')- (''') 



In the equation above N(ijj) = {^^^I^bT _ j^^-i jg ^.j^g 
erage number of reservoir thermal excitations , with fcs 
the Boltzmann constant and T the reservoir tempera- 
ture, and J(w) is the spectral density of the environment 
defined, in the continuum limit, as 



J(^) 



(8) 



with m„ the masses of the environmental oscillators. In 
deriving the master equation no Markovian approxima- 
tion has been done. The memory effects are included in 
the time-dependent coefficients A(<) and 7(t). The latter 
term is known as dissipation coefficient and gives rise to 
a classical damping term that is not dependent on tem- 
perature. The former term A(t) is known as diffusion 
coefficient and is directly proportional to the reservoir 
temperature [1]. 

It is worth mentioning here that performing a secular 
approximation does not affect the non-Markovian short 
time dynamics of certain observables in the weak cou- 
pling limit |l6j. In this paper we focus on the dynamics 
of one of such observables, namely the heating function. 
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III. MODELING DIFFERENT TYPES OF 
RESERVOIRS 

We now introduce a class of spectral densities in or- 
der to compare the QBM dynamics for different types of 
reservoirs. The spectral densities we examine are of the 
form 



(9) 



The exponential cutoff is introduced to eliminate diver- 
gencies in the a; — > oo limit. We have compared various 
types of cutoff functions and concluded that their dif- 
ferent analytical forms do not play a major role in the 
dynamics of QBM. For the sake of simplicity we therefore 
focus, in the rest of the paper, on the exponential cutoff. 
The parameter s appearing in Eq. © is a constant that 
can acquire values < 1, 1 or > 1, corresponding to the 
so called sub-Ohmic, Ohmic and super-Ohmic spectral 
densities, respectively. In this paper we consider some 
examples and fix the value of s to 1/2, 1 and 3. The 
three cases have different physical interpretations. The 
Ohmic spectrum [s — 1) gives, for QBM, a friction-like 
force that is proportional to velocity. The Ohmic spec- 
trum can be used, e.g., to describe charged interstitials 
(conductive electrons) in metals [2I]. The super-Ohmic 
spectral density (s = 3) corresponds to, e.g., a phonon 
bath in one or three dimensions, depending on the sym- 
metry properties of the strain field ^ . It is also possible 
to show that this type of environment can be used in de- 
scribing the effect of the interaction between a charged 
particle and its own electromagnetic field [1^. The sub- 
Ohmic spectral density (s = 1/2 ) corresponds to the 
type of noise that may occur in some solid state devices 
and, in the high T case, is similar to the "1/f noise" in 
Josephson junctions [i^ . 

We introduce the spectral distribution given by 



/(w) = J{uj) 



N{u) + i 



(10) 



This quantity contains all the information needed, in the 
weak coupling limit, about the reservoir, i.e., the density 
of modes and the occupancy of each mode. The spectral 
distribution depends on the temperature of the reservoir 
through the average number of reservoir thermal excita- 
tions N(uj). At high temperatures T we can approximate 
N{u}) « ksT /uj while at zero temperature N(u)) = 0. In 
the rest of the paper we discuss the QBM dynamics in 
these two temperature regimes. 

A relevant parameter in the description of QBM is the 
resonance parameter r defined as the ratio between the 
cutoff frequency ujc and the frequency of the system os- 
cillator Wo, i.e.. 



Wo' 



(11) 



We consider three exemplary values of the resonance pa- 
rameter, namely r equal to 0.1, 1 and 10. The r = 0.1 



case is characterized by the fact that the effective cou- 
pling between the system oscillator and the environment 
is very small for all three reservoir types because the 
system oscillator is detuned from the peak of the reser- 
voir spectral distribution. We call this the off-resonant 
case. The Ohmic and sub-Ohmic reservoirs are such that 
the effective coupling between the system and the reser- 
voir becomes stronger when r grows from 0.1 to 10. The 
super-Ohmic reservoir, on the contrary, shows the highest 
effective coupling for r = 1, while the r = 0.1 and r = 10 
cases correspond to relatively weak couplings. Plots of 
the spectral distribution for different values of r in the 
high and zero temperature limits are shown in Figs. [1] 
andm 

The spectral distribution at high temperatures, given 

by 



/(w) = a^keT 



/u 



(12) 



is shown in Fig. [T] The sub-Ohmic spectrum has a diver- 
gency point at w = 0. This causes large effective coupling 
induced by the low frequency part of the spectral density. 
An opposite example can be found for the super-Ohmic 
reservoir, where the peaks of the spectrum in the cases 
r — 1 and r — 10, lie in the higher frequency range. 
The spectral density for zero T reservoirs, given by 



J(w) 



(13) 



is shown in Fig. [5] (the plots are here grouped according 
to parameter r for clarity) . A key difference with respect 
to the high temperature case is that at zero tempera- 
ture, the sub-Ohmic spectrum does not diverge in zero 
anymore. 

As we will see in Sees. IV and V, the QBM heat- 
ing dynamics depends crucially on the different form of 
the spectral distributions in the Ohmic, sub-Ohmic and 
super-Ohmic cases. Our main goal is to establish a clear 
connection between the reservoir properties and the dy- 
namics of both the decay rate and the heating function, 
as done, for an Ohmic reservoir with Lorentzian cutoff, 
in Ref. ji^]. In this way we will be able to motivate 
from a physical point of view the origin of the different 
QBM dynamics for different reservoirs, and therefore in 
different physical contexts. 



IV. DECAY RATES 

The frontfactors [A(i) -|- 7(t)]/2 and [A(i) -7(t)]/2 in 
the master equation (O represent the relaxation rates for 
the two decay channels of the QBM model. In the Fock 
state basis the former rate is associated to \n) \n — 1) 
transitions, i.e. to the transfer of one excitation from the 
system to the environment. The latter rate corresponds 
to \n) — > |rt + 1) transitions, i.e. it describes the ab- 
sorption of one excitation from the environment. These 
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FIG. 1. Spectral distributions of different reservoirs at high 
temperatures. Here / — I / {o?kBT) and u = lo/luo- For each 
spectral curve the location of the cutoff frequency is given by 
tJc = i^c/oJo = The location of the oscillator frequency has 
been marked with a solid vertical line. 



transitions, describing heating or cooling of the quantum 
harmonic oscillator due to the interaction with the ex- 
ternal environment, destroy the quantum coherence of 
initial superpositions. 

After a certain reservoir-dependent time, the decay 
rates reach their constant positive Markovian values 



Am = 7r/(i:jo), 
7m = ■^-'"('^o)- 



(14) 
(15) 



The decay rates can also temporarily attain negative 
values. When this happens, the corresponding decay 
channel has been shown to operate in a reverse way. 
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FIG. 2. Spectral distributions at zero temperature grouped 
according to the resonance parameter r. Here I — 2//(a^tJo) 
and Lj = lj/luq. The vertical dotted-dashed line is the location 
of the cutoff frequency ijJc = cJc/i^o- 



i.e., the down channel actually induces heating and vice 
versa [45| . We will now examine the decay rates in non- 
Markovian time scales at high and zero temperatures. 



A. High temperatures 

At high temperatures A(t) ^ 7(^)1 ^-^d for time scales 
much shorter than the thermalization time, both the 
transition up and down channels operate at approxi- 
mately the same rate A{t)/2. We have obtained an ana- 
lytic expression for A{t) = A(t)/ {aksT) for all the three 
types of environment considered in this paper. More pre- 
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FIG. 3. Decay rates at high temperatures for different types of reservoir in the non-Markovian time scales. Here A 
A/(2a^fcsr) and r = 0.1 
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FIG. 4. The decay rates at zero temperature for r — 10. The figure on the left depicts the transition up rates, while the figure 
on the right corresponds to the transition down rates. 



cisely, for the Ohmic environment we obtain 



environment is 



'i cosh 



z + 2t 



sinh 



SI (^) - SI 



r 

z + 2t 
r 



(16) Mt) 



and for the super-Ohmic environment 



At cos I 



2 sin I 



(l+t2)2 r + rt^ 
't 



(17) 



i cosh 



ci 



z + 2t 



sinh 



St (^) - St 



z + 2t 



where z = i — t, ci{x) and si(x) are the cosine and sine 
integrals defined as ci{x) = — ^^^^^^dx and si{x) — 
— /J^ ^^^^^dx. Finally, the decay rate for the sub-Ohmic 



lr{l+it) 



X + erf(z") - ie^/'Vl + erf(iz") (18) 

+ i^Jl^I^t [erf (z+) - ie2/''erf(iz+) 



[(1 + i)^fi±l\l y/2^ and erf(a;) = 



From the analytic expressions of the diffusion coeffi- 
cient A(t) one can show that such a quantity oscillates 
taking temporarily negative values for r 1 for all reser- 
voir types. The origin of these oscillations can be traced 
back to the form of the spectral distribution. Figure 1 
shows that, e.g., for r = 0.1 most part of the spectrum 
is located in a region of the frequency space such that 
a; < wq. Figure 3 shows the oscillatory behavior of A(t) 
for this value of r. As r grows, the spectrum starts to 
overlap with wq. We note that, in the case of the sub- 
Ohmic environment, the diffusion coefficient presents os- 
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cillations for all values of r, but it attains negative values 
only when r <^ 1. 

The Markovian value of the diffusion coefficient A(<) is 
proportional to I{luq) [See Eq. p3) ]. As one can see from 
Fig. [l] for r = 10 the highest value of Am corresponds 
to the sub-Ohniic environment while in the off-resonant 
r = 0.1 case Am is small for all reservoir types. 



B. Zero temperature 

When the system oscillator interacts with a zero- 
temperature reservoir, the Markovian theory predicts 
that energy is transferred from the system to the envi- 
ronment, i.e., the oscillator is driven towards its ground 
state. Our non-Markovian theory, however, shows that 
for times ujct ^ 1, the average energy of the system oscil- 
lator may increase, as we will see in detail in Sec. V. This 
is due to the form of the interaction Hamiltonian, given 
by Eq. ([3]) , containing four terms characterizing the emis- 
sion and absorption processes, namely, a6„, a&Jj, a^6„ 
and a^6jj. The two terms in the middle correspond to 
real processes conserving the unperturbed energy, while 
the other two are known as the counter rotating terms. 
These terms describe the simultaneous creation or anni- 
hilation of a quantum of energy both in the system and 
in the reservoir oscillators. The energy required for such 
processes to occur comes from the system-reservoir cou- 
pling. By combining these two counter rotating terms, 
we obtain a process that corresponds to an energy con- 
serving process. It has been shown that at zero temper- 
ature the dynamics of the decay rate for the transitions 
up originates from these counter rotating terms [l7j . 

The decay rates at zero temperature show similar de- 
pendence on the parameter r as the one discussed in the 
high temperature case. In particular, for r <C 1 both 
A{t) + 7(t) and A{t) — ^{t) oscillate attaining negative 
values, a clear signature of the non-Markovian behavior 
of the system. Compared to the high T case, oscillations 
in A{t) — 7(<) obtaining negative values are present also 
for higher values of r, e.g. r = 1, for all reservoir types. 
For the super-Ohmic environment the decay rates show 
a strong initial jolt for all values of r. Having in mind 
the form of the spectral distributions one sees that the 
initial jolt is present whenever the peak of the spectrum 
lies in the frequency region u) > lOc- 

For ujct » 1 the decay rate A{t) — j{t), describing 
\n) — > |n + 1) transitions in the system oscillator, ap- 
proaches zero as expected from the Markovian theory, 
while A{t) + 7(t), describing |n + 1) — > \n) transitions, 
reaches a constant positive value Am +"/m = 7r/(ajo), as 
shown in Fig. |4]for the r = 10 case. 

An interesting aspect worth mentioning is visible in the 
super-Ohmic case for r = 10. Figure |3] clearly shows that 
for the super-Ohmic spectrum the decay rate A{t) +j{t), 
i.e. the cooling rate of the system oscillator, approaches 
its small but nonzero constant value already for ujct ~ 3, 
while at the same time the decay rate A{t) — 7(i), de- 
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FIG. 5. Thermalization times for different reservoirs. Here 
I 2 , 

tth = TTQ tth- 



scribing heating of the system oscillator, attains negative 
value and tends to zero while remaining negative. 

It has been shown in Ref. (451] that in correspondence 
to negative regions of the time dependent coefficients, 
in our case A(t) ±"f(t), reverse transitions restoring the 
previous quantum state occur. In view of these results 
one can argue that, in the case considered above, the 
up channel acts like a transition down channel. This 
implies that the thermalization is achieved via a reverted 
transition up channel, while the actual transition down 
channel is almost completely closed. 

In the next section we will see how the behavior of 
the decay rates is related to the heating dynamics of the 
quantum Brownian particle and we will investigate the 
differences in the dynamics of the mean energy of the 
system due to different environments. 



V. HEATING OF A QUANTUM BROWNIAN 
PARTICLE 



A. Markovian thermalization dynamics 

The QBM dissipative dynamics can be described by 
means of the heating function, defined as 



{n) = a'a. 



(19) 



The analytical expression for the heating function is given 

by Mm 



{n{t)) ^ e~^'^'\n{0)) + 



1 r 



-m 



1 



+ Ar(i), (20) 



where r{t) and Ar{t) are defined as 



T{t) = 2 / -f{h)dh 
Jo 

ft 



Ar{t) = e-r« / e^^*'^Ait,)dh 
Jo 



(21) 
(22) 
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In the following we focus on the case where the os- 
cillator is initially in the ground state, i.e., (n(0)) — 0, 
and perform a comparative study of the heating function 
dynamics for different reservoir structures. 

We begin by looking at the Markovian dynamics de- 
scribing the time evolution for times much greater than 
the reservoir correlation time. The Markovian expression 
of the heating function {n(t))M is obtained inserting the 
Markovian expressions of the diffusion and dissipation 
coefficients, given by Eqs. (HI])- (HI]), into Eqs. (pg)) -(P ^ . 

(n(i))M-A^(c^o)(l-e""), (23) 

where T — 2jm = 7rJ(wo)- From Eqs. © and (ITSI) we 
can express the reservoir thermalization time, in units of 
Wo, as follows 

= ujo/T = i7ra^)-\'-'e'/\ (24) 

From Fig. [5] we see that, for both the sub-Ohmic and 
the Ohmic reservoirs, the thermalization time decreases 
monotonically for increasing values of r, i.e. for increas- 
ing values of the cutoff frequency with respect to the 
frequency of the system oscillator ujq. On the contrary, 
for the super-Ohmic reservoir, there exist a value of r, 
namely r « 0.5, minimizing the thermalization time. In 
general, all the three reservoir types considered in this pa- 
per are such that the thermalization time grows rapidly 
when r — >■ 0, and correspondingly the thermalization pro- 
cess is notably slowed down. Our analysis suggests that, 
by appropriately changing the cutoff frequency of a high 
temperature engineered reservoir, it is possible to control 
the thermalization dynamics. 

B. Non- Markovian heating 

We now look at the non-Markovian short time dynam- 
ics of the heating function. For times much smaller than 
the thermalization time Eq. (|20p can be approximated 
by 

{n{t)) = / [A(ti)-7(ti)] dti. (25) 
Jo 

This equation establishes a clear connection between the 
heating function dynamics and the time dependent decay 
rate A(i) — j{t) corresponding to transitions increasing 
the system oscillator energy, i.e., absorption of quanta 
from the environment. This is clearly related to our 
choice of the initial condition {n{0)) = 0. In this case 
indeed, for times much shorter than the thermalization 
time, the absorption of a quantum of energy from the 
environment (heating) dominates over the opposite pro- 
cess, i.e., the emission of a quantum of energy into the 
environment (cooling). 



1. High Temperatures 

For high temperature reservoirs, Eq. ((25|) can be fur- 
ther approximated by (47j 

(nit)) « / A{t,)dt,. (26) 
Jo 

The sign of A{t) determines whether the heating function 
grows monotonically or exhibits an oscillatory behavior. 
In more detail, when A{t) oscillates taking negative val- 
ues then the heating function oscillates. 

The oscillations in the heating function appear when 
the system gives back to the reservoir some of the en- 
ergy that had previously been absorbed from it. In other 
words, the direction of the energy flow is reversed during 
the time periods in which the slope of {n{t)) is negative. 
These oscillations are a sign of the non-Markovian dy- 
namics and are due to the finite reservoir memory that 
allows the partial and temporary recovery of some of the 
information/energy lost in the reservoir. 

Equation ((26|) links the heating function to the diffu- 
sion coefficient dynamics. In the previous section we have 
seen how the structure of the reservoir, and in particu- 
lar some specific system-reservoir parameters, determine 
the temporal behavior of A{t). By means of Eq. ((26)) 
we can now establish a connection between the reser- 
voir spectrum and the heating process and compare the 
non-Markovian heating for the sub-Ohmic, Ohmic and 
super-Ohmic reservoirs. 

Prom Fig. 6 one sees clearly that the non-Markovian 
heating, described by Eq. ()26p . presents two main types 
of behavior: an oscillatory behavior and a monotonic 
growth similar to the Markovian heating. In Ref. [i^ 
it was demonstrated that, for an Ohmic reservoir and for 
r = 0.1, oscillations in {n{t)) originate from the low fre- 
quency part of the spectrum while the monotonic heating 
is caused by the resonant part of the spectrum, namely 
by the value of /(wq). 

This connection between the features of {n{t)) and the 
spectrum seems to hold also for the three reservoir types 
discussed in this paper. Indeed, for r > 1 the decay rates 
A{t) are always positive giving rise to monotonic heating. 
Moreover the bigger is the value of I{ujo), the faster is the 
non-Markovian monotonic heating rate. For r = 1, e.g., 
the monotonic heating occurs at approximately the same 
rate for all the reservoir types corresponding to the fact 
that /(ojo) is the same for all the reservoirs. 

The connection between the low frequencies of the 
spectrum and oscillations in the heating dynamics can 
be illustrated by considering the cases where oscillations 
are present, i.e., for r <^ 1. This is the parameter re- 
gion where A{t) obtains temporarily negative values giv- 
ing rise to oscillations in {n{t)). For these values of r 
the peaks of all the three spectral distributions are posi- 
tioned in the low frequencies region, indicating that the 
presence of oscillations and the low frequency part of the 
spectrum are intertwined. 
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FIG. 6. Short time dynamics of the heating function for different high temperature reservoirs. In the plots a = 0.01 and we 



have set ^ = 100. 



Oscillations in {n{t)) mark the presence of non- 
Markovian memory effects. The persistence of non- 
Markovian effects for a given value of r (e.g., r <C 1 for 
high T reservoirs) depends on the type of reservoir spec- 
trum. Figure 6 shows that the memory effects persist 
for much longer times in the sub-Ohmic reservoir than 
in the other reservoir types. This indicates that, when 
dealing with sub-Ohmic reservoirs, a non-Markovian ap- 
proach might be needed also at time scales where, for 
the Ohmic and super-Ohmic environments, a Markovian 
treatment is sufficient. 

We have now illustrated the main features of the non- 
Markovian heating dynamics of QBM for high tempera- 
tures. In the next section we will briefly summarize the 
corresponding results for the zero temperature case. 



efficient A{t) — 7(i), associated to the heating rate of the 
system oscillator, tends to zero for times much greater 
than the reservoir correlation time, as expected from the 
Markovian theory. After the initial non-Markovian heat- 
ing, the system will eventually thermalize with the zero 
temperature reservoir. 

Similarly to the case of high T reservoirs, also at low 
temperatures oscillations in {n{t)) vanish as r becomes S> 
1. Now, however, this typically non-Markovian feature is 
present for greater values of r than in the high T case, 
e.g., also for r — 1. The fact that non-Markovianity 
is present for a wider range of r, in the case of a zero 
temperature reservoir, is in agreement with what found 
m Ref. fi^ for a Ohmic reservoir with Lorentz-Drude 
cutoff. 



2. Low Temperatures 

When the system oscillator, initially prepared in its 
ground state, interacts with a zero temperature reser- 
voir, the dynamics of the heating function is basically 
due to the finite, although small, system-reservoir cou- 
pling energy. At zero temperatures, and for short times, 
the heating function dynamics is linked to the decay co- 
efficient A{t) - 7(i) via Eq. pS]). 

The non-Markovian dissipative dynamics of QBM at 
zero T is qualitatively similar to the one shown in Fig. 6 
for high T reservoirs. As discussed in Sec. IIVBI the co- 



VI. CONCLUDING REMARKS AND FUTURE 
PROSPECTS 

In this paper we have presented analytical results on 
the dynamics of a quantum Brownian particle interacting 
with different types of bosonic thermal reservoirs. Our 
approach makes it possible to identify in a clear way the 
microscopic physical processes taking place at short non- 
Markovian times scales and to link them to the spectral 
properties of the environment. By comparing the effects 
that different types of environments have on the quantum 
system, we can identify those physical contexts where a 
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non-Markovian approach is required in the description of 
the time evolution. Moreover we can predict the type 
of reservoir that perturbs less the quantumness of the 
system of reference, in our case the harmonic oscillator, 
and the value of the parameters for which this occurs. 
We have seen indeed that the decay and heating rates 
for the QBM depend strongly on the reservoir type, and 
in particular on the ratio r between the cutoff frequency 
LOc and the frequency of the system oscillator ujq. 

The time dependent coefficients appearing in the mas- 
ter equation and in the heating function are connected via 
Eqs. (PS)) and (1^51) for both zero and high T reservoirs. 
The heating dynamics for both temperature regimes 
shows a similar dependence on the parameter r. In gen- 
eral oscillations in the heating function, a typical non- 
Markovian feature, are present for r <C 1. For zero tem- 
peratures, however, these oscillations are also present for 
r « 1. So for all the types of environment considered, 
zero T reservoirs are inherently more non-Markovian 
than high T reservoirs. Moreover, we demonstrated that 
the sub-Ohmic reservoir induces more pronounced and 
longer lasting non-Markovian dynamics compared to the 
other reservoir types and gives rise to a faster Markovian 
heating in the resonant case (r = 10). 

The oscillations in (n(t)) indicate a back and forth ex- 
change of energy between the system and the reservoir. 
In particular, when the the slope of {n{t)) is negative, 
the system gives back some of the excitation received 
previously from the reservoir. The statistical ensemble 
of system oscillators is, in general, colder than the reser- 
voir with whom it interacts, since we assumed through- 
out the paper that the system is prepared in its ground 
state. Oscillations in the heating function therefore in- 
dicate the possibility of a temporary energy flow from a 
cooler object to a hotter one due to the memory of the 
environment. 

This is not a surprise. The monotonic evolution of the 
system density operator to its equilibrium value, which is 
a universal property of quantum dynamical semigroups 
(Spohn's theorem) is in general violated at short 

(non-Markovian) timescales. This anomaly has been pro- 
posed to be used to control the thermodynamics of an 
atomic system simply by changing the way in which it is 
measured 49]. 

Our results indicate that by means of reservoir engi- 
neering techniques, e.g., by changing the parameter r, 
one could modify the thermalization dynamics of the sys- 
tem. Another intriguing possibility stemming from reser- 
voir engineering is the simulation of paradigmatic models 
of open quantum systems as the one discussed in this pa- 
per. Schemes for simulating QBM with trapped ions were 
presented in Refs. 22, 42]. A similar approach may be 
used to simulate the sub-Ohmic and super-Ohmic envi- 
ronments here investigated. 

An ideal physical context where our results could be 
experimentally verified is the trapped ion context. Ex- 
periments with single trapped ions have demonstrated 
the ability to engineer artificial environments and to 



control the relevant system-environment parameters [6|. 
These experiments aim at measuring the decoherence of a 
quantum superposition of coherent states and Fock states 
due to the presence of the reservoir. Several types of en- 
gineered reservoirs are demonstrated, e.g., thermal am- 
plitude reservoirs, phase reservoirs, and zero temperature 
reservoirs Q. 

A high T amplitude reservoir is obtained by apply- 
ing a random electric field E whose spectrum is centered 
on the axial frequency lOz of oscillation of the ion. The 
trapped ion motion couples to this field due to the net 
charge q of the ion: Hint — ~qx ■ E, with x — {X, Y, Z) 
displacement of the cm. of the ion from its equilibrium 
position. Remembering that E oc ei{bi + bl), with hi 

and bl annihilation and creation operators of the fluctu- 
ating field modes, and that X oc (^a + a^) the quantized 
position operator of the ion motion, one realizes that this 
coupling is equivalent to the bilinear one given in Eq. ([3]) . 

The random electric field is applied to the endcap elec- 
trodes through a network of properly arranged low pass 
filters limiting the "natural" environmental noise but al- 
lowing deliberately large applied fields to be effective. 
This type of drive simulates an infinite-bandwidth am- 
plitude reservoir It is worth stressing that, for the 
times of duration of the experiment the heating due to 
the natural reservoir is definitively negligible @ . 

The different high T spectra that we discuss in this 
paper can be realized experimentally by filtering the ran- 
dom field, used in the experiments for simulating an 
infinite-bandwidth reservoir, with a suitable set of band- 
pass filters. This enables the comparison of the heating 
rates between different reservoir spectra and the observa- 
tion of non-Markovian effects. We notice that measure- 
ments of the heating function are routinely performed in 
the trapped ion context. The heating function is obtained 
from measurements of the population of the vibrational 
states of the ion. A detailed study of the experimen- 
tal techniques to simulate harmonic quantum Brownian 
motion with trapped ions, for the case of an Ohmic spec- 
trum with Lorentz-Drude cutoff has been performed in 
Ref. [l^l . We believe the methods needed to verify these 
phenomena are already in the grasp of the experimental- 
ists [Ulil]. 

Finally, it is worth recalling that manipulation of the 
dynamics via certain types of measurements causes the 
decay processes to be inhibited or accelerated, depend- 
ing on the system-reservoir properties. These crucially 
quantum phenomena, known as quantum Zeno (QZE) 
and anti-Zeno effects (AZE), have been studied for the 
QBM model in the case of an Ohmic spectrum in Ref. 

The results presented in this paper pave the way to the 
study of the infiuence of the reservoir spectrum on the 
occurrence of quantum Zeno or anti-Zeno effect (AZE). 
The borderline between the occurrence of QZE or AZE is 
indeed related to the spectral properties of the environ- 
ment and therefore will depend on the type of reservoir 
(Ohmic, sub-Ohmic, super Ohmic) considered. A com- 



10 



parativc study of the Zeno-anti-Zeno crossover is thus a 
natural follow up of this paper. 
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